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1—5 Abstract 

pg We significaiitly advance the research program initiated in "Moment-Based Evidence for Sim- 

1""! pie Rational- Valued Hilbert-Schmidt Generic 2x2 Separability Probabilities" (J. Phys. A, 45, 

'j^ 095305 [2012]). A function P{a), incorporating a family of six hypergeometric functions, all with 

H argument g| = (|)^, is obtained here. It reproduces a series, a = g, 1, |, 2, . . . ,32, of sixty- 

^-( four Hilbert-Schmidt generic 2x2 generalized separability probabilities, advanced on the basis of 

systematic high-accuracy probability-distribution-reconstruction computations, employing 7,501 

(N 

^ determinantal moments of partially transposed 4x4 density matrices. For generic (9-dimensional) 

^ two-rebit systems, P{-^) = g|, (15-dimensional) two-qubit, P{1) = ^, and (27-dimensional) two- 



quat(ernionic)bit systems, P{2) = oJr. The function P{a) is generated-applying the Mathematica 



cn 

O command FindSequenceFunction-usmg solely eUher the thuty-two integral or half-mtegral proba- 



bilities, yet successfully predicts the unused thirty-two. 
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I. INTRODUCTION 

The predecessor paper [T] -addressing the relatively long-standing 2x2 separability 
probability question [2HII] (cf. [T^HTi] )-consisted largely of two sets of analyses. (In 
the foundational study [2], "three main reasons of importance" -philosophical, practical 
and physical-were given for studying this question.) The first set of analyses in |1] was 
concerned with establishing formulas for the bivariate determinantal product moments 
P^^\ \p\ ) ,k,n = 0,1,2,3, ... , with respect to Hilbert-Schmidt (Euclidean/flat) mea- 
sure [T3] [ini sec. 14.3], of generic (9-dimensional) two-rebit and (15-dimensional) two-qubit 
density matrices (p). Here p^^ denotes the partial transpose of the 4x4 density matrix 
p. Nonnegativity of the determinant Ip"^"^] is both a necessary and sufficient condition for 
separability in this 2x2 setting [T7] . 

In the second set of analyses in [1], the univariate determinantal moments (Ip^^l ) 
and ((|p^"^| IpD"), induced using the bivariate formulas, served as input to a Legendre- 
polynomial-based probability distribution reconstruction algorithm of Provost [ TBI sec. 2]. 
This yielded estimates of the desired separability probabilities. (The reconstructed proba- 
bility distributions based on \p^'^\ are defined over the interval |p^"^| G [—jq, ^], while the 
associated separability probabilities are the cumulative probabilities of these distributions 
over the nonnegative subinterval \p^'^\ € [0, i^]. We note that for the fully mixed (classical) 
state, |p^"^| = 2565 while for a maximally entangled state, such as a Bell state, \p^'^\ = —jq-) 

A highly-intriguing aspect of the (not yet rigorously established) determinantal moment 
formulas obtained (by C Dunkl) in [1, App.D.4] was that both the two-rebit and two-qubit 
cases could be encompassed by a single formula, with a Dyson-index-like parameter a [19] 
serving to distinguish the two cases. The value a = \ corresponded to the two-rebit case 
and a = 1 to the two-qubit case. (Let us note that the results of the formula for a = 2 
and n = 1 and 2 have recently been confirmed computationally by Dunkl using the " Moore 
determinant" (quasideterminant) [201 EI] of 4 x 4 quaternionic density matrices. However, 
tentative efforts of ours to verify the a = 4 [conjecturally, octonionic [22j, problematical] 
case, have not proved successful.) 

When the probability-distribution-reconstruction algorithm [18] was applied in [1] to the 
two-rebit case (a = ^), employing the first 3,310 moments of Ip''^^], a (lower-bound) estimate 
that was 0.999955 times as large as || ~ 0.453120 was obtained (cf. [231 p. 6]). Analogously, 



in the two-qubit case {a = 1), using 2,415 moments, an estimate that was 0.999997066 
times as large as Jj ~ 0.242424 was derived. This constitutes an appeahngly simple rational 
value that had previously been conjectured in a quite different (non-moment-based) form of 
analysis, in which "separability functions" had been the main tool employed [ll| . 

Further, the determinantal moment formulas advanced in pj were then applied with a 
set equal to 2. This appears-as the indicated recent (Moore determinant) computations of 
Dunkl show-to correspond to the generic 27-dimensional set of quaternionic density matrices 
[2HI2S]- Quite remarkably, a separability probability estimate, based on 2,325 moments, that 
was 0.999999987 times as large as j^ ~ 0.0804954 was found. (In line with this set of three 
results, the paper [Ij was entitled, "Moment-Based Evidence for Simple Rational- Valued 
Hilbert- Schmidt Generic 2x2 Separability Probabilities".) 

In the present study, we extend these three (individually-conducted) moment-based anal- 
yses in a more systematic, thorough manner, jointly embracing the sixty-four integral and 
half- integral values a = ^, 1, |, 2, . . . , 32. We do this by accelerating, for our specific pur- 
poses, the Mathematica probability-distribution-reconstruction program of Provost [18j, in 
a number of ways. Most significantly, we make use of the three-term recurrence relations for 
the Legendre polynomials. Doing so obviates the need to compute each successive higher- 
degree Legendre polynomial ab initio. 

In this manner, we were able to obtain-using exact computer arithmetic throughout- 
" generalized" separability probability estimates based on 7,501 moments for a = 
|, 1, |, . . . , 32. In Fig. Q^we plot the logarithms of the resultant sixty-four separability prob- 
ability estimates (cf. P Fig. 8]), which fall close to the line -0.9464181889a. In Fig. [2] we 
show the residuals from this linear fit. 

Most notably, in Fig.[3]we will present a hypergeometric-function-based formula, together 
with striking supporting evidence for it, that appears to succeed in uncovering the functional 
relation (P(a)) underlying these generalized separability probabilities. Further, in (pi), and 
the immediately preceding text, we list a number of remarkable values yielded by the hyper- 
geometric formula for values of a other than the basic 64 (half-integral and integral) values 
from which we have started. 
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FIG. 1: Logarithms of generalized separability probability estimates, based on 7,501 Hilbert- 
Schmidt moments of |p^^|, as a function of the Dyson-index- like parameter a 
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FIG. 2: Residuals from linear fit to logarithms of generalized separability probability estimates 



II. RESULTS 



A. The three basic conjectures revisited 



1. a = 2~the two-rebit case 



In [T] , a lower-bound estimate of the two-rebit separability probability was obtained, with 



the use of the first 3,310 moments of \p \. It was 0.999955 times as large as 



29 
64 



0.453120. 



With the use, now, of 7,501 moments, the figure increases to 0.999989567. This outcome, 
thus, fortifies our previous conjecture. 



2. a = 1-the two-qubit case 

In [T], a lower-bound estimate of the two-qubit separability probability was obtained, 
with the use of the first 2,415 moments of Ip'^^l, that was 0.999997066 times as large as 
|j ^ 0.242424. Employing 7,501 moments, this figure increases to 0.99999986. 

3. a = 2~the quaternionic case 



In [1], a lower-bound estimate of the quaternionic separability probability was ob- 

26 
323 



tained that was 0.999999987 times as large as t^^ ?« 0.0804954, using the first 2,325 mo- 



ments of Ip^^I- Based on 7,501 moments, this figure increases, quite remarkably still, to 
0.999999999936. 

B. Generalized separability probability formula 

A principal motivation in undertaking the analyses reported here-in addition, to further 
scrutinizing the three specific conjectures reported in [T]-was to uncover the functional 
relation underlying the curve in Fig. [l] (and/or its original non- logarithmic counterpart). 

Preliminarily, let us note that the zeroth-oidei approximation (being independent of the 
particular value of a) provided by the Provost probability-distribution-reconstruction al- 
gorithm is simply the uniform distribution over the interval [—jq, ^]- The corresponding 
zeroth-order separability probability estimate is the cumulative probability of this distribu- 
tion over the nonnegative subinterval [0, ^], that is, ^/(^^ + ^) = 17 ~ 0.0588235. So, 



it certainly appears that speedier convergence (sec. II A ) of the algorithm occurs for separa- 
bility probabilities, the true values of which are initially close to ^ (such as t^ ~ 0.0804954 
in the quaternionic case). Convergence also markedly increases as a increases. 

It appeared, numerically, that the generalized separability probabilities for integral and 
half-integral values of a were rational values (not only for the three specific values a = 
|, 1, 2 of original focus). With various computational tools and search strategies based upon 
emerging mathematical properties, we were able to advance additional, seemingly plausible 
conjectures as to the exact values for a = 3,4,..., 32, as well. (We inserted many of our 
high-precision numerical estimates into the search box on the Wolfram Alpha website-which 
then indicated likely candidates for corresponding rational values.) 



We fed this sequence of thirty-two conjectured rational numbers into the FindSequence- 
Function command of Mathematica. (This command "attempts to find a simple function 
that yields the sequence a^ when given successive integer arguments.") To our considerable 
satisfaction, this produced a generating formula (incorporating a diversity of hypergeomet- 
ric functions of the pFp_i type, p = 7, . . . , 11, all with argument z = || = (|)^) for the 
sequence. (Let us note that 2;" 2 = . /|| is the "residual entropy for square ice" [261 P- 412] 
(cf. pTj eqs.[(27), (28)]) [28]. In fact, the Mathematica command succeeds using only the 
first twenty-eight conjectured rational numbers, but no fewer.) 

However, the formula produced was quite cumbersome in nature (extending over several 
pages of output). With its use, nevertheless, we were able to convincingly generate rational 
values for /ia//-integral a (including the two-rebit || conjecture), also fitting our correspond- 
ing half-integral thirty-two numerical estimates exceedingly well. (Let us strongly emphasize 
that the hypergeometric-based formula was generated using only the integral values of a. 
The process was fully reversible, and we could first employ the half-integral results to gen- 
erate the formula-which then-seemingly perfectly fitted the integral values.) 

At this point, for illustrative purposes, let us list the first ten half-integral and ten in- 
tegral rational values (generalized separability probabilities), along with their approximate 
numerical values. 



a = i I 0.453125 " = 1 i 0.242424 

a = I ^^ 0.137562 a = 2 ^ 0.0804954 

" = i liiiii 0.0480083 « = 3 liH 0.0290081 

« = i JIMll^ 0.0176948 « = 4 Jfi, 0.0108722 

a = I i°5i504368i347 0.00671852 a = 5 J^m^ 0.00417101 



2 9007199254740992 "-""^'-^^"^ -^ " 21460999 

1 

a 



11 r9256°294^^^^^^^^^^^^^ 0.0025994 a = 6 .}I!IT^?.^ 0.00162519 



2 27670116110564327424 "-""^"^^^ '^ " 110638410169 

^13 3387374833367307236269 fl 001 fll SQ n = 7 191151001 () r\r\r\f,Amr\Q 

" 2 3324546003940230230441984 u.uu±u±o» u ( 298529164591 u.uuuu^iuouy 

^15 124792688228667229196729 n nnnAn'iOOJ ry = R 1331199762 n nnnOKA'^q^ 

2 309485009821345068724781056 u.uuu'±uoz,z, i u o 5232880523393 u.uuuz,o^u»j. 

n- = 11 407557367133399293946182513 000160753 rt- = 9 74195568677 000101799 

2 2535301200456458802993406410752 "•">-"-' -l^'J ' ^^ " ^ 729345064647247 '^•^'-"-'-l^-l ' ^^ 



, = 19 1338799759394288468677657208071 o 0000644609 ry = 10 730710456538 o 00004089' 

2 20769187434139310514121985316880384 U.UUUUU'±'±UU3 w lu 17868447453498669 ^-^^^^^^^^^ 

(1) 

To simplify the cumbersome (several-page) output yielded by the Mathematica FindSe- 



quenceFunction command, we employed certain of the " contiguous rules" for hypergeometric 
functions listed by C Krattenthaler in his package HYP [21] (cf. [SO])- Multiple applications 
of the rules C14 and C18 there, together with certain gamma function simplifications sug- 
gested by C Dunkl, led to the rather more compact formula displayed in Fig. pi (Attempts 
to achieve a still more succinct form have not yet succeeded.) This formula incorporates 
a six- member family {k = 1, . . . , 6) of hypergeometric functions, differing only in the first 
upper index k, 

^ / 2 3 4 5 7 6 13 3 17 19 21 27 

rFgl fc,a +-,« + -,« +-,« + -,« + -,«+-;« + —,« + -,«+—,«+— ,a + 2,a+— ; — 

(2) 
Interestingly, we are only able to, in general, evaluate the formula numerically, but then to 
arbitrarily high (hundreds, if not thousand-digit) precision, giving us strong confidence in 
the validity of the exact generalized separability probabilities that we advance. Whether 
exact values can be directly extracted, without recourse to such high-precision numerics, 
appears to be an issue still yet to be fully resolved. 

Let us now apply the formula (Fig. [3]) to values of a other than the basic sixty-four. For 
a = 0, the formula yields-as would be expected-the "classical separability probability" of 1. 
Further, proceeding in a purely formal manner (since there appears to be no corresponding 
genuine probability distribution over [—jq, ^]), for the negative value a = —^, the formula 

yields |. For a = — ^, it gives -2. Remarkably still, for a = |, the result is clearly (to 

i7r(i)^ 
one thousand decimal places) equal to 2 "^ — ^v = 2 — „, rk'^i/-, ~ 0.6486993992, where 

^ ' ^ 21agm(l,v^) 2l^/2■n■i/^ ' 

the arithmetic-geometric mean of 1 and -\/2 is indicated. (The reciprocal of this mean is 
Gauss's constant.) For a = f , the result equals 2 — ,,„„ r-J't5\ ~ 0.3279684732, while for 

a = -I, we have ^^^^^ + 2 = 2+ ^^T^^J-* ^ 7.087249321. For a = |, the outcome is 

4' 21agm(l,v^) 2l7r'3/2 3' 
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FIG. 3: Hypergeometric formula P{a) for Hilbert-Schmidt generic 2x2 generalized separability 
probabilities and evidence that it reproduces the basic three (real [a = 5], complex [a = 1] and 
quaternionic [a = 2]) conjectures of ||, Jj and ^ 



288927r(i) 
3440807r2 



0.36424897456. Results are presented in the table 



a 

_3 

4 

_2 
3 

_1 
2 

_1 
3 

_1 

4 

1 

4 

1 
3 

2 
3 



2 + 
2- 



2 + 

2- 
2- 



P(a) 

32V2r(|^ 

2l7r3/2' 
Btt 

2 
3 

3rrif 



2 

2lV27r3/2 
459v^7r 

2 _ 288927r(i) 



3 
\ 4 



2- 



3440807r2 

9689r(|) 

4420VSr(|) 



value \ 

7.08725 
1.24527 

0.666667 

3.461 

2 

0.648699 
0.572443 

0.364249 
0.327968 



(3) 



(Let us note that the term 



3r 



47r2 



1.46099848 present in the result for a 



■| is "Baxter's 



four-coloring constant" for a triangular lattice j26| p. 413].) Also, for a = —1, we have |. 



For a 



■|, the result is 



2 
'3' 



Thus, it appears that our hypergeometric-related formula (Fig. pi) constitutes a major 
advance in addressing the long-standing (Hilbert-Schmidt) 2x2 separability probability 
question [2j . However, there certainly remain the important problems of formally verifying 
this formula (as well as the underlying determinantal moment formulas in p], employed in 
the probability-distribution reconstruction process), and achieving a better understanding 
of what it conveys regarding the geometry of quantum states (TU]. Further, questions of 
the asymptotic behavior of the formula {a — )• oo), possible modified/simplified/canonical 
forms of it, and of possible Bures metric [31 HI El [161 El] counterparts, are under current 
investigation. 
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